Armchair graphene nanoribbons: Electronic structure and electric field modulation 
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We report electronic structure and electric field modulation calculations in the width direction 
for armchair graphene nanoribbons (acGNRs) using a semi-empirical extended Hiickel theory. Im- 
portant band structure parameters are computed, e.g. effectives masses, velocities and bandgaps. 
For the three types of acGNRs, the p z orbital tight-binding parameters are extracted if feasible. 
Furthermore, the effect of electric field in the width direction on acGNRs dispersion is explored. It 
is shown that for the two types of semiconducting acGNRs, an external electric field can reduce the 
bandgap to a few meV with different quantitative behavior. 

PACS numbers: 73.22.-f, 73.20.-r, 72.80.Rj 



I. INTRODUCTION 

Unconstrained graphene is a two dimensional hexag- 
onal monolayer of carbon atoms. Its unique linear dis- 
persion around the Dirac point and zero bandgap^ has 
generated significant interest 3 -^. Constraining one di- 
mension of graphene results into nanoribbons. The elec- 
tronic structure of these graphene nanoribbons (GNR) 
depends on the width and chirality 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ■ Two 
unique GNRs are armchair and zigzag referred to as 
acGNR and zzGNR in this article. acGNR has an arm- 
chair edge as shown in Fig. 1 and when conceptually 
rolled to form a nanotube results in a zigzag tube and 
vice versa. Some experimental techniques have already 
been used to measure their properties^ 3 - and numerous 
fabrication schemes have been devise d 15 i 16 ' 17 . Electronic 
applications of graphene and GNRs are also being sought 
afteri^^. 

In zzGNRs, the wavefunctions for conduction and va- 
lence bands are localized at the edges^. In addition, 
the bands around the Fermi energy have very small dis- 
persion that leads to Stoner magnetism^. These edge 
states can be modulated with an external electric field in 




FIG. 1: Electronic structure of armchair graphene nanorib- 
bons (acGNR). The ball and stick model of a graphene 
nanoribbon with N=9 is shown with the unit cell. E-k di- 
agrams are shown for three different types of acGNRs using 
extended Hiickel theory (EHT). 



the width direction, resulting in half mctallicity 21 . 

In acGNRs, the wavefunctions associated with bands 
around Fermi energy are distributed throughout the 
width of the nanoribbon. However, these bands still 
can be modulated with an external electric field in the 
width direction as discussed by Novikov using a contin- 
uum model^. In addition, due to quantization in one 
direction, acGNRs have velocities less than those found 
in unconstrained graphene sheets, and the band structure 
has a parabolic character around the band edge within a 
few tens of meV. 

In this paper, we focus on acGNRs and study their 
electronic structure and electric field modulation in the 
width direction with a semi-empirical extended Hiickel 
theory (EHT). Similar electric field modulation effects 
have been studied in carbon nanotubes as wel l 23 ' 24 . The 
detailed model has been reported in Refi^. EHT pa- 
rameters are transferable and have been benchmarked 
with generalized gradient approximation of density func- 
tional theory (DFT) for carbon atoms in graphene struc- 
ture. EHT is computationally inexpensive and hence 
appropriate for calculating properties of large systems 
without compromising accuracy. As an example, up 
to about 1000-atom electronic structure calculations 2 ^ 
and up to about 150-atom transport calculations 2 ^ have 
been reported in silicon based systems with modest com- 
putational resources. In this paper, up to about 160 
atoms calculations are presented. Contributions from 
five nearest neighbors are included. C-C atomic dis- 
tance is taken as 1.44A, for which EHT parameters have 
been optimized. We find that incorporating about 3.5% 
decrease in C-C atomic distance^ 7 , for the edge carbon 
atoms results in a bandgap increase of about 52meV for 
N=18 W=2.1nm and bandgap decrease of about 64meV 
for N=19 W=2.2nm acGNR. Small bandgap changes 
are expected since wavefunctions for valence and con- 
duction bands are deloaclized for acGNR. However, in 
zzGNR, these wavefunctions arc localized on the edges 
and any atomic relaxation would have significant effect 
on the bandstructure. Since these variations are small 
in acGNR, we ignore any atomic relaxation in the re- 
ported calculations. Atomic visualization is done using 
GaussVicw2£. 
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FIG. 2: (color online) Bandgaps and effective masses, (a) Variation of bandgap with nanoribbon widths of different types of 
acGNRs. Using a p z -orbital tight binding method, t=2.5eV and t=2.7eV match the bandgaps obtained by extended Huckel 
theory (EHT) for /3-acGNRs and 7-acGNRs, respectively, (b) Variation of effective mass with nanoribbon widths of different 
types of acGNRs. 



II. ELECTRONIC STRUCTURE 

On a p z level of the tight-binding theory, two thirds 
of acGNRs are semiconducting with a bandgap in- 
versely proportional to their widths and the other has 
zero bandgap depending on the chirality^. However, 
one obtains a different result using the more sophisti- 
cated theory^, such as EHT and DFT. First, the zero 
bandgap acGNRs also have a small bandgap that is in- 
versely proportional to the width. Second, the remain- 
ing semicondcuting acGNRs only follow an inverse re- 
lation within its own category. For convenience, we 
propose to categorize them into a-, (3- and 7-acGNRs. 
This classification is similar to the ones used recently 
in Rcfs^i^. a-acGNRs are N=8,ll,14,... and have 
very small bandgap. /3-acGNRs are N=9,12,15,... and 
7 acGNRs are N=10, 13,16,... acGNRs have also been 
classified into three subclasses in context of the orbital 
diamagnetismi. 

An electronic structure calculation for each type of 
acGNR is shown in Fig. 1. As can be seen that 
N=8 a-acGNR has a small bandgap and has a nonlin- 
ear dispersion around the T point. N=9 /3-acGNR has a 
large bandgap with a parabolic dispersion around the T 
point. Interestingly, N=10 7-acGNR has a slightly larger 
bandgap with larger effective mass dispersion around the 
r point and smaller velocity in the linear region away 
from the the T point as compared to N=9 /3-acGNR. 
We extract the bandgaps and effective masses within a 
few tens of meV around the band edges of these three 
types of acGNRs and plot them in Figs. 2(a) and (b), 
respectively. Fig. 2(a) is a computational verification of 
earlier results^ on a semi-empirical level. We find that 
incremental change in the bandgap of 7-acGNRs with 
respect to /3-acGNRs is smaller in EHT than local den- 
sity approximation of density functional theory^. For 
each type of acGNR, bandgaps and effective masses are 
inversely proportional to the width with a different pro- 
portionality constant. The bandgap versus width (W) 



relations are given as: 

r QMeV/W{nm) for a - acGNR 
E gap = \ 0.86eV/W(nm) for (3 - acGNR 
{ 1.04eV/W(nm) for 7 - acGNR 

We find Fig. 2(b) important because some ap- 
proaches toward graphene structures involve effective 
mass description^. Each type of acGNRs follow an in- 
verse relation of effective mass with the width given be- 
low: 

r 0.005/W(nm) for a - acGNR 
— = \ 0.091/W(tito) for f3 - acGNR 
m ° { 0.160/W{nm) for 7 - acGNR 

where m is the free electron mass. It should be noted 
that using a p z -orbital tight binding model, the effective 
mass follows the same inverse relation versus width for 
all three types of acGNRs^. Furthermore, we determine 
the p z -orbital tight binding parameters that reproduce 
the bandgaps as shown in Fig. 2(a). These parameters 
are 2.5eV and 2.7eV for j3- and 7- acGNRs respectively. 
Since tight-binding parameter for 7-acGNRs is higher, 
we conclude that wavefunctions are hybridized more in 
this type of acGNR. This physical effect has some impli- 
cations for electric field modulation as discussed in next 
section. 



III. ELECTRIC FIELD MODULATION 

Fig. 3(a) shows electric field modulation of the band 
structure for an N=10 7-acGNR. The effective mass 
around the T-point increases with increasing electric field 
(E) and eventually changes sign, similar to Refi^. Fur- 
thermore, for E=0, the band dispersion in the linear 
regime away from the T-point shows velocity very close 
to the unconstrained graphene velocity (=8.8 x 10 5 m/s) 
indicated by red (grey) circles. With increasing E, veloc- 
ity in this linear regime away from the T-point decreases 




FIG. 3: (color online) Electric field modulation of band dispersions, (a) Variation of velocity in the width direction for N=10 
7-acGNR. The linear dispersion shown by red circles represents a value of 8.8 x 10 m/s - velocity around the Dirac point 
for graphene calculated using EHT. (b) Variation of effective masses. Effective masses are obtained by parabolic fits to the 
conduction bands within a few ksT of band edge for /3-acGNRs and 7-acGNRs, and within a fraction of a ksT for a-acGNRs. 
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FIG. 4: Bandgap modulation. Bandgap as a function of width and electric field for (a) /3-acGNRs and (b) 7-acGNRs. 7-acGNRs 
have larger bandgap modulation as compared to /3-acGNRs. 



to about 5 x 10 5 m/s. In addition, the bandwidths of 
the valence and conduction bands are also decreasing. 
Moreover, a Mexican hat structure is observable that has 
been seen in acGNR^ 2 -, carbon nanotubes 2 ^ and graphene 
bilayer o 9 ' 30 ! 31 ' 32 . These features are in qualitative agree- 
ment with the electric field effects reported in semicon- 
ducting acGNRs elsewhere using a continuum model^ 2 -. 
However, there are some quantiative differences which we 
address in this section. We show the extracted effective 
masses around the T-point for N=8, 9 and 10, which are 
a-, (3- and 7-acGNRs, respectively in Fig. 3(b). These 
effective masses arc valid for tenths of ksT for a-acGNRs 
and for a few UbT for (3- and 7-acGNRs. After this en- 
ergy scale, the band dispersions become linear again and 
remain so for about a few electron volts when they be- 
come nonlinear and hence saturate as shown in Fig. 3(a). 
lV/nm electric field is within the dielectric breakdown 
limit of thermal Si02 , which may result in higher electric 
field inside graphene due to smaller dielectric constant. 
Moreover, high-K dielectrics can be used to further en- 
hance the electric field. However, such a high electric 
field may lead to dielectric reliability issues and is unde- 



sirable. 

In addition, the bandgap is modulated with increasing 
electric field. A clear feature is the location of wavevec- 
tor corresponding to the conduction/ valence band mini- 
mum/maximum. These two perturbations in the band 
structure arc further explored in Figs. 4 and 5 re- 
spectively for (3- and 7-acGNRs. In Fig. 4, we show 
bandgap modulation as a function of width and electric 
field. A threshold behavior is observed, similar to Refi^ 2 -, 
where bandgap starts decreasing appreciably above a 
threshold electric field E t . The dimensionless parame- 
ter u t = eE t x W/Eg a p is reported as 4.5 for both kinds 
of acGNRs in Refi^ using a continuum model. How- 
ever, we find that this is different for these two acGNRs 
and is about 5.6 and 3.9 for (3- and 7-acGNRs respec- 
tively. Moreover, for 7-acGNRs, bandgap decreases at a 
faster rate compared to /3-acGNRs and thus 7-acGNRs 
have larger bandgap modulation. This is consistent be- 
cause wavefunctions are more hybridized in 7-acGNRs 
and hence any perturbation affects the band structure 
more than /3-acGNRs. 

Additionally, the bandgap for /3-acGNRs monotoni- 
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FIG. 5: (color online) The wavevector corresponding to band 
field for (a) /3-acGNRs and (b) 7-acGNRs. The value of k at 
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cally decrease with electric field. However, the bandgap 
decreases appreciably only after a threshold electric field. 
This is different from Refi^ 2 -, where below the threshold 
electric field, bandgap is constant and it decreses only 
after the threshold electric field. Furthermore, for 7- 
acGNRs, the bandgap first increases a little and then 
decreases - a feature although small, but not present in 
continuum calculations 22 . 

With an appropriate electric field applied, one can re- 
duce the bandgap of a semi-conducting acGNR to a few 
mcV. We find that bandgap never becomes zero, whereas 
using a continuum models, one finds zero bandgap. In 
order to change the band structure, one has to incorpo- 
rate perturbation on the order of the tight-binding pa- 
rameter (2.5eV for (3- and 2.7eV for 7-acGNRs). There- 
fore, an electric field of lV/nm should not be able to 
induce a significant change in small width acGNRs due 
to small perturbation as shown in Fig. 4. However, the 
same electric field can change the electronic structure of 
a wider acGNRs due to larger potential variation. The 
physics behind this bandgap narrowing is the spectral 
shift of the conduction and valence band states on the 
two edges. This leads to downward and upward shift 
for conduction and valence band, respectively. Further- 
more, in Fig. 5, we show the wavevector shift (ksE) cor- 
responding to conduction band minimum/ valence band 
maximum. Again, 7-acGNRs have larger shift as com- 



pared to /3-acGNRs. Overall, this shift can be as much as 
one third of the wavevector at X point. Unfortunately, we 
could not find a consistent set of p z -orbital tight-binding 
parameters to reproduce Figs. 3, 4 and 5 simultaneously. 



IV. CONCLUSIONS 

We have studied band structure and electric field 
modulation of acGNRs using EHT. The three types of 
acGNRs exhibit distinct electronic structure and electric 
field modulation properties. We extract important band 
structure parameters and a set of p^-orbital tight-binding 
parameters benchmarked with extended Hiickcl theory to 
reproduce the bandgaps. Additionally, electric field mod- 
ulation results are compared with a continuum model 2 ^. 
We find that qualitative trends are the same, however 
there are some quantitative differences between the two 
models. 
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